Non-adiabatic dynamics in 10 Be with the microscopic a+a+iV+iV model 
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The a+ 6 He low -energy reactions and the structural changes of 10 Be in the microscopic 
a+a+N+N model are studied by the generalized two-center cluster model with the Kohn-Hulthen- 
Kato variation method. It is found that, in the inelastic scattering to the a+ 6 He(2 1 ') channel, char- 
acteristic enhancements are expected as the results of the parity-dependent non-adiabatic dynamics. 
In the positive parity state, the enhancement originates from the excited eigenstate generated by the 
radial excitation of the relative motion between two a-cores. On the other hand, the enhancement 
in the negative parity state is induced by the Landau-Zener level-crossing. These non-adiabatic 
processes are discussed in connection to the formation of the inversion doublet in the compound 
system of 10 Be. 
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In the last two decades, developments of experiments 
with secondary RI beam have extensively proceeded the 
studies on light neutron-rich nuclei. In particular, much 
efforts have been devoted to the investigation of molec- 
ular structure in Be isotopes. Theoretically, molecular 
models with 7r and a orbitals along the axis connecting 
two a-particles have been successful in understanding the 
low- lying states of this isotopes Q, ■ Experimentally, 
the molecular structures were mainly investigated by the 
breakup processes 0, and the sequential decays [|| us- 
ing the high energy RI beams. 

In recent experiments, furthermore, the low-energy 
6 He beam becomes available. Low-energy reaction cross- 
sections such as the elastic scattering with an a target 
0, Q and the sub-barrier fusions with heavy target [3 
have been accumulated. In future experiments, it will 
also be possible to investigate the molecular states in Be 
isotopes through the reactions such as a+ 6 He and a+ 8 He 
with low-energy 6,8 He-beams. Therefore, it is very inter- 
esting to study theoretically on the low-energy scattering 
of 6 He and 8 He by an a target. 

In studying reaction processes exciting the molecular 
degrees of freedom, it is very important to construct a 
unified model which is capable of describing both struc- 
ture and reaction on the same footing. For this purpose, 
we introduce a microscopic model, the generalized two- 
center cluster model (GTCM) 0, El In this model > ^ 
is possible to describe both molecular and atomic limit 
of the system of C1+C2+N+N+... where (7, is the i-th 
cluster core and N is the nucleon. In the region where 
two core nuclei are close, the total system is expected to 
form the molecular orbital structure, while in the region 
where two core nuclei are far apart, the molecular orbitals 
smoothly change into product wave functions consisting 
of the atomic orbitals. 

In this paper, we apply the GTCM for the 10 Be nucleus 
with the a+a+N+N four-body model. We will ana- 
lyze both molecular structure in 10 Be and the low-energy 
a+ 6 He scattering. Besides the description of the a+ 6 He 
reaction, such analysis will be useful to understand the 
breakup mechanism of 10 Be into clusters. In spite of 



many theoretical efforts in the last decade P,0jE3 _ E3> 
only Ref. [HI discusses the molecular-orbital formation 
in 10 Be and the a+ 6 He scattering problem in a unified 
way. 

Current experimental investigations are extended to 
the 12 Be and 14 Be nuclei 0, U and hence, theoretical 
studies extended to such heavier systems are urged. The 
GTCM approach has a potential to 12 Be and 14 Be as 
well as 10 Be. On the other hand, the direct extension 
of the approach in Ref. becomes quite difficult for 
systems with many valence neutrons, since it utilizes the 
a+a+N+N few-body model. 

One of the favorable features of our approach is that 
it is possible to describe the nuclear Landau-Zener (L-Z) 
transition microscopically ^f|-[2(j. The L-Z transition 
is induced by the avoided level-crossing of two adiabatic 
potential-surfaces. When the avoided crossing occurs in 
the potential surfaces, the adiabatic states drastically 
change their intrinsic character at the crossing point. For 
large relative velocity, the colliding nuclei follows not the 
adiabatic path but the diabatic one. This diabatic mo- 
tion was discussed by Landau and Zener for the atomic 
collision [2l|. The nuclear analogue has been called the 
nuclear L-Z transition. In spite of long history on this is- 
sue, a clear evidence of the L-Z transition is still lacking 
in nuclear collision [2^, ■ Our microscopic approach 
will indicate a possible nuclear L-Z transition in the col- 
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First, we briefly explain the framework of GTCM 0, 
IToj . The basis functions for 10 Be are given as 
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A{ip L (a)^ R (a)tp(m)ip(n)} . (1) 



The ev-cluster wave function ipi( a ) (i—L,R) is given by 
the (0s) 4 configuration in the harmonic oscillator (HO) 
potential. The position of an a-cluster is explicitly spec- 
ified as the left (L) or right (R) side. The relative motion 
between a particles is described by a localized Gaussian 
function specified by the distance S [24[ . A single-particle 
state for valence neutrons around one of a clusters is 
given by an atomic orbitals (AO), tp(i,pi-,T) with the 
subscripts of a center i (-L or R), a direction p/- (k=x, 



y, z) of Op-orbitals and a neutron spin r (=f or J.). In 
Eq. (JIJ, the index m(n) is an abbreviation of the AO 
(i,Pk,r). The intrinsic basis functions with the full anti- 
symmctrization A are projected to the eigenstate of the 
total spin J, its intrinsic angular projection K and the 
total parity tt by the projection operator . 

The total wave function is finally given by taking a 
superposition over S and K as 

^ = fdS CRS) <S>f K {S) (2) 

with i = (m, n). The coefficients Cf(S) are determined 
by solving a coupled channel GCM (Generator Coordi- 
nate Method) equation [24| . If we fix the generator coor- 
dinate S and diagonalize the Hamiltonian with respect to 
i and K, we obtain the energy eigenvalues as a function 
of S, which we call the adiabatic energy surfaces (AES). 

In the present calculation, we used the Volkov No. 2 and 
the G3RS for the central and the spin-orbit part of the 
nucleon-nucleon (NN) interaction, respectively. The pa- 
rameters in the NN interactions are modified from those 
in Ref. ;9] so as to reproduce the threshold of a+ 6 He 9 . s . 
and the excitation energy of the 6 He(2+) state . This 
is because the reproduction of the threshold is essentially 
important in the treatment of the scattering problem. 

The Majorana (M), Bartlett (B) and Heisenberg (H) 
exchanges in the central part are fixed to M=0.643, 
B=— 11=0.125, while the strength of the spin-orbit force 
is chosen to 3000 MeV for the repulsive part and 2000 
MeV for the attractive part. The radius parameter b of 
HO wave functions for a clusters and valence neutrons 
is commonly taken as 1.46 fm. We included all the AO 
configurations of two neutrons that can be constructed 
by the Op-orbitals. 

The AES with J n =0+ is shown in Fig. □ There ap- 
pear three local minima at the short distance region of 
the AES. The adiabatic states (AS) at the lowest and 
third minima have the molecular orbital configurations 
of {it~i 2 ) 2 and (7r^ 2 ) 2 ' respectively The AS 

around the second minimum has a dominat configura- 
tion Of (cr+ 2 ) 2 nm, but the one particle excited con- 
figuration, (&f, 2 Trf, 2 ), * s a ^ so stron gly mixed. The latter 
configuration has a spin triplet structure and hence, the 
coupling between them is induced by the two-body spin- 
orbit interaction. 

It is well known that the simple crjy 2 orbital is not suf- 
ficient to describe both the Oj" state in 10 Be and the l/2 + 
state in 9 Be as discussed in Ref. lj. Itagaki et al. shows 
that spin-orbit interaction generates the strong coupling 
betweeen the distorted (c^) 2 configuration with the 

spin-triplet configuration and the pure (c^ 2 ) 2 one, which 

plays very important role for lowering the 0^" state in 
10 Be pj. Therefore, the present result is consistent to 
that discussed in Ref. pj and it is reasonable to describe 
the intrinsic structure of 10 Be(0 2 f ). 
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FIG. 1: Adiabatic energy surfaces for J n —0 + . The low- 
est and third local minima have a dominant configurations of 
(V3/2) 2 and (tt-^) 2 , respectively, while the second one has a 
configutaion of the distorted (o"^/ 2 ) 2 configuration (See text 
for details). The lowest and the second surfaces with the 
double circles are the dinuclear states of [a+ 6 He(0^)]L=o and 
[a+ 6 He(2] l ")]_L = 2, respectively, while those with the solid cir- 
cles have the configurations of [ 5 He(3/2~)+ 5 He(3/2~)] with 
(iX) = (00) and (22), respectively. 
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FIG. 2: The same as Fig.Qbut for the negative parity states 
(J w — 1~). The dashed (solid) surfaces has a dominant compo- 
nent of {^/2 a t/2) K=1 (( 7r i/2 cr i/2) jr = 1 ) aroun d the local mini- 
mum, while it is smoothly connected to the dinuclear channel 
of [a+ 6 He(2^")] with L — 1 [L = 3) at an asymptotic region. 
The surface with a solid circles has a dominant component of 
the [a+ 6 Re(0f)] L =i channel. 



At the asymptotic region (S > 8 fm) where two a- 
cores are completely separated, the valence neutrons are 
localized around one of the a cores. The localization 
of the orbitals leads to the formation of the dinuclear 
channels such as [ 4 He+ 6 He(/)]L (double circles) and 
[ 5 Hc(/i)+ 5 Hc(/ 2 )]/l (solid circles), in which individual 
channels are specified by the intrinsic spin of the clusters 
(I 1; I 2 ), the channel spin / (I=Ii + 12) and the relative 
angular momentum between clusters, L. The asymptotic 
energy position of the lowest AES is higher by about 5 
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FIG. 3: Energy spectra for J 7r =0 + . The error bar means the 
decay width of the resonance states. The solid curve shown 
in the right part is the partial cross section for the inelastic 
scattering [a+ 6 He s . 3 .]L=o — > [a+ 6 He(2^)]_L = 2 with a respec- 
tive scale at the up-most axis. The AES in Fig. are shown 
by the dashed curves. 



MeV than the cv+ 6 He g s . threshold. This is because the 
relative motion between clusters in Eq. (|TJ is described 
by the locally peaked Gaussian |24| and hence, its kinetic 
energy contributes to the AES in the asymptotic region. 

The structural changes occur smoothly between the 
molecular orbital region and the dinuclear channels re- 
gion in passing through the intermediate region shown by 
the arrow in Fig. ^ jf| ^j. In the intermediate coupling 
region, we can see the level crossing between the sur- 
face of (^3/2) 2 (dashed curve) and the second one (c^) 2 
(solid curve). The energy splitting at the crossing point 
is about 1.5 MeV. 

In Fig. [3 the AES for the J n =l~ state is shown. The 
configurations of the valence neutrons smoothly changes 
in the AES for the a— a distance parameter except for the 
curves with the solid circles. The AS along this surface 
has an almost pure-component of the [a-r- 6 He(0^)],L = i 
channel. Thus, this AES is not molecular orbitals but 
the dinuclear state in the whole regions of a — a distance. 
Because of the different character of the lowest two or- 
bitals, the AES behaves different as a function of a — a 
distance. In contrast to the results of the J 71 " = + states 
shown in Fig. ^ this causes an clear avoided-crossing at 
S=6 fm as shown by the dotted circles in Fig. [21 The 
energy splitting at the crossing point is about 0.5 MeV 
which is smaller than that in J 7r =0 + (~1.5 MeV). This 
means that, in the negative parity state, the change of 
the intrinsic structure is much sharper in the distance 
than the case of the positive parity state. 

To take into account the excitation of the relative mo- 
tions between two a-cores, we solve the GCM equation 
by employing the AS from S=l fm to S=70 fm with 
the mesh of 0.5 fm. The calculated energy spectra of 
jtt = q+ an( j jtt = i~ are shown in Figs. |3| and respec- 
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FIG. 4: The same as Fig. but for J*=l~. In the right 
part, the partial cross section for the inelastic scattering 
of [a+ 6 He 9 . s .]L=o — > [Q+ 6 He(2 1 ")]L=i is shown by the solid 
curve. 



tively. In solving GCM, we apply the Absorbing-kernels 
in the Generator Coordinate Method (AGCM) in which 
the absorbing boundary condition is imposed outer re- 
gion of the total system [2j|. Due to the absorbing 
boundary, the resonance poles can be clearly identified 
in the complex energy plane. 

In the positive parity, the 0^~, 0^ and 0|" states are the 
poles corresponding to the respective local minima in the 
AES, having the molecular orbital structures. Therefore, 
we should call these states the "adiabatic poles" , because 
they can be realized as the local minima in the AES. On 
the other hand, there is no local minimum corresponding 
to the 0^" state and it is generated by the radial excitation 
on the distance parameter S. Thus, the 0^" state should 
be called as the "radially-excited poles", although it is 
the pole generated by a linear combination of the AS. 
The structure of Oj" is very different from other lower 
adiabatic poles, because it is orthogonalized to the adia- 
batic poles and is spatially extended. The wave function 
in 0^ has an enhanced component of [a+ 6 He(2 ] f )]l=2 
at the surface region and hence, it is different from the 
molecular orbital configuration. 



In J w =l , we have identified two adiabatic poles cor- 
responding to the two local minima of i^J^i/^K^x- The 
lower pole (J z — 3/2) is generated from the linear combi- 
nation of the AS around the lowest local minimum, while 
the higher one (J z = 1/2) is originated from the minimum 
in the third AES as indicated by the arrow in Fig. 0] 

In the present calculation, we identify no resonance 
corresponding to the Q+ 6 He s . s . AES, although its ap- 
pearance is discussed within the bound state approxima- 
tion in Ref. @. Therefore, the a+ 6 He £ ,. s . cluster config- 
uration will not be stabilized as a resonance pole under 
the present condition reproducing the respective thresh- 
old energy. Since the a+ 6 He g . s . AES has a flat shape 
over the wide range of the distance Q, the stablility of 
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this configuration becomes quite sensitive to the energy 
position of the threshold. 

Let us consider the difference in the J v —0 + and 1" 
in relation to the inversion doublet structure. The 
(ir 3 



. 3/2 ) 2 AES in J*=0+ and the a- 



- 6 Re g . s . AES in J n =V 



originally form the inversion doublet of a+ 6 He s . s . with 
A' 7r =0 ± . The J 7r =0 + partner couples to the symmetric 
5 He+ 5 He configuration and hence, it is strongly distorted 
to the molecular orbital as the distance S gets closer. In 
the J n =l~ surface, however, the dinuclear configuration 
of a+ 6 He g s is well developed, which is similar situation 
to the 20 Ne=a+ 16 O system \2($. Therefore, the appear- 
ance of the avoided crossing in the negative parity has a 
close connection to the formation of the inversion doublet 
in the compound system of 10 Be. 

We next discuss how the AES profile and the level 
scheme appear in the low energy reactions. We show 
the partial cross section for the inelastic scattering to 
6 He(2^) excitation. In solving the scattering prob- 
lem, we employ the Kohn-Hulthen-Kato (KHK) vari- 
ation method [23 where the AO basis in Eq. 
are transformed to the asymptotic channel wave func- 
tion. The present calculation is equivalent to the usual 
coupled-channels calculation including the channels of 
a+ 6 He(0^,0^,2+,2+,l+) and 5 He(3/2-,l/2-)+ 5 He(3/ 
2~,l/2~). Our calculation thus includes much more 
channel components than the previous studies of a+ 6 He 
cluster model |12(. In this calculation, the maximum S is 
changed from 70 fm to 12 fm and the channel wave func- 
tion is matched to the scattering Coulomb wave function 
at a matching radius fm. 

The calculated partial inelastic cross-sections are 
shown in the right part of Figs. and In J w =0 + , 
a strong peak appears at E c , m ~ 3 MeV, although there 
is no definite avoided crossing in the AES. This is due to 
the effect of the radially-excited pole, 0^", which include 
the large component of the exit 2^ channel. We can also 
see the enhancements at E c . m , ~ 7 MeV, which nicely 
coincident to the adiabatic pole of (tt^) 2 ' O4 , but it is 
much smaller than that generated by the radially-excited 
poles. 

In J n =l~ , the strong enhancement can also be seen at 
Ec.m.^ 6 MeV, nevertheless there is no pole in the inci- 
dent and exit channels. The energy of the enhancement 
is quite close to that of the avoided crossing at S=6 fm. 
To investigate the origin of this enhancement, we have 
solved the coupled-channel problem between the lowest 
two AS by employing the adiabatic Kohn-Hulten-Kato 
method that will be explained in the next paragraph. 
In such calculation, we have confirmed that this peak is 
generated by the coupling between the lowest two AS, 
nevertheless there appears no poles. Therefore, we can 
conclude that this inelastic peak is due to the L-Z tran- 
sition at the avoided crossing. Though the (^jy 2 CT i/2) 
resonance is located close to the cross section peak, it 
weakly couples to the incident and exit channels. This 
pole is found to just generate the kink at a slightly lower 
energy than the peak position. 
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FIG. 5: The partial cross section of the a+ 6 He 9 . s . elastic 
scattering. The left and right panels show the results for 
J w —0 + and J w — 1~, respectively. The dotted curve shows 
the result of the adiabatic approximation, while the solid and 
dashed ones show that of two AS coupled-channel and that 
of the full coupled-channels, respectively, 



To see the connection between the AS and the scat- 
tering process in a transparent way, we formulate the 
adiabatic KHK (AKHK) method in which individual AS 
are employed as the basis functions in solving the scat- 
tering problem. In the following, we briefly explain the 
formulation of the AKHK. First, we define the a-th AS 
at a distance S by 



iK 



> J cl(S) 



x X l(R,S)} ((3 = hI 2 IL) , 



(3) 



(4) 



where ^^ nK (S) is the AO basis given by Eq. @J. In 
the second line, the AS is expanded in terms of the 
channel (CH) wave function, ^^^(S). Eq. |@} shows 

the explicit expression of ^ch(^) which is constructed 
from the angular momentum coupling among the inter- 
nal states of i-th nucleus Vij.Afj an d the spherical har- 
monics Y^jw(R) with the relative coordinate R. In the 
last line of Eq. J3J, Xl(R, S) denotes the locally peaked 
Gaussian with the peak position of R ~ S. 

The mixing coefficients in the a-th AS, F^(S), satisfies 
the following relation 



hm F$(S) 



I0,a 



(5) 



Eq. JSJ means that the a-th AS becomes a specific chan- 
nel a at an asymptotic region (S — > 00), although a vari- 
ous channel components are strongly mixed in the inter- 
nal region (S ~ small). In solving the scattering problem 
with the basis of the AS, therefore, only a specific chan- 
nel a satisfying Eq. JSJ) should be transformed into the 
scattering basis function, because the respective radial 
function xl {R, S) do not satisfy the scattering boundary 
condition. 
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The transformation can be done by utilizing the 
KHK method, in which the localized basis functions 
are smoothly connected to the scattering Coulomb wave 
function at a matching radius Rc [I?! • The Rc should be 
taken to be sufficiently large value where all the channel 
components except for a are completely damped. The 
calculation of the matrix element for the AKHK basis 
can be easily done by the similar procedures shown in 
Ref. [27|. The details of the AKHK method will be shown 
in a forthcoming paper. 

In order to discuss the gross features of the non- 
adiabatic effects, we calculate the partial cross sections of 
the elastic scattering. In two panels of Fig. [5J the results 
calculated from the AKHK method are shown. In this 
calculation, the S m ax and the Rc are takne to be the 
same values in calculating the inelastic cross section. In 
both panels, the dotted curves show the pure adiabatic 
approximation. That is, only the AS along to the lowest 
AES is employed in solving the scattering problem. The 
solid curves show the solution of the full coupled channel 
(CC) in which all the AS are included. 

In the result of J 7r =0+, the adiabatic approximation 
quite nicely simulates the full CC solution in the low- 
energy region. Furthermore, the adiabatic approximation 
simulates the gross behavior of the full CC solution up to 
about 4 MeV except for the kink just below 2 MeV. This 
means that the elastic scattering mainly proceeds along 
to the lowest AES. The dashed curve shows the result in 
which the lowest two AS are coupled. The coupling effect 
improves the adiabatic approximation. The dashed curve 
deviates from the solid one in the higher energy than 3 
MeV, but the difference is not so large. 

The results of J 7r = l _ is drastically different from those 
of J n =0 + . The validity of the adiabatic approximation is 
limited only in the region of the small cross section. Fur- 
thermore, the non-adiabatic coupling between the lowest 
two AS strongly reduces the cross section of the adiabatic 
approximation, which amounts to about one order re- 
duction. Therefore, the adiabatic approximation is quite 
poor for describing the scattering process. The main con- 
tribution of the non-adiabatic coupling is from the first 
excited AS, because the dashed curve is similar to the 



solid one. This is due to the appearance of the avoided 
crossing which can be clearly seen in Fig. [3 

In summary, we investigated the adiabatic proper- 
ties of the 10 Be—a+a+N+N structures as well as 
the a+ 6 He g . s . low-energy reactions, especially the non- 
adiabatic dynamics including the Landau-Zener (L-Z) 
transition. We achieved such an unified study of the 
structures and the reactions in the generalized two- 
center cluster model (GTCM). The scattering problem 
in the GTCM is solved by applying the Kohn-Hulthen- 
Kato (KHK) variation method and the adiabatic KHK 
(AKHK) method proposed in the present study. We have 
found the non-adiabatic enhancements in the inelastic 
scattering cross-section to the a-|- 6 He(2^) in both the 
JT—Q+ anc j states. However, their origins are 

very different to each other. In the J"*=Q + , the enhance- 
ments is due to the appearance of the radially-excited 
pole in the 10 Be system. This is realized by the radial 
excitation of the a-core's relative motions. In contrast, 
the enhancement in J 7T = l~ is originated from the L-Z 
level-crossing between the lowest two adiabatic energy 
surfaces. Such a difference of the non-adiabatic dynamics 
also affects the elastic scattering process. In the positive 
parity, the adiabatic approximation is good for describing 
the collision process, while it becomes wrong approxima- 
tion in the negative parity. 

Next important step is to compare our result with 
recent experiments 0, by extending the present ap- 
plication to the higher partial waves. In a comparison 
with experiments, we should be careful to optimize the 
micleon-nucleon interaction so as to reproduce the ob- 
served energy spectra in 10 Be together with the thresh- 
old of a+ 6 He and 5 He+ 5 He. Such extended studies are 
under progress. 
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